Dynamics of small perturbations against a stationary polariton distribution 
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The dynamics of small perturbations against the stationary density matrix of a pumped polariton 
system with only one photon polarization is studied. Depending on the way the system is pumped 
and probed, decay times ranging from 30 to 5000 ps are found. The large decay times under resonant 
pumping are related to a bottleneck effect in the decay of the excess (probe) populations of dark 
polariton states. No singular behaviour at the threshold for polariton lasing is observed. 
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Polariton lasers arc lasers without population inver- 
sion p~] . Coherence buildup in them is the result of the 
quasibosonic statistics of the polaritons, i.e. quasipar- 
ticles composed by electron-hole pairs from a quantum 
well strongly interacting with photons from a semicon- 
ductor microcavity. They share similarities with ordi- 
nary photon lasers and Bose-Einstein condensates [21 [3] . 
In these devices, the threshold power for lasing is 1 - 2 
orders of magnitude lower than in ordinary lasers. Room- 
temperature polariton lasing has been recently reported 
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Besides these promising characteristics, polaritons in 
microcavities provide an exceptional possibility for fun- 
damental research. An example is the recent paper [S], 
where the authors seek for evidences of the Goldstone bo- 
son, which should appear in connection with the buildup 
of coherence in the polariton system. Indeed, in a cylin- 
drical microcavity, the two (approximately) circular po- 
larizations of the fundamental photon mode are related 
to two degenerate polariton "condensates" . The relative 
phase between the two photon polarizations acts as an 
order parameter. Phase fixing leads to a linear polariza- 
tion [5J, whose direction can be easily rotated (the Gold- 
stone excitation). In paper 5 , Ballarini et. al. study 
the changes in the PL response of a planar microcav- 
ity induced by small pulse perturbations of the pumping. 
They measure the lifetime of these perturbations, show- 
ing that it is much higher than the cavity decay time, 
and that it increases when the stationary pumping ap- 
proaches the threshold for polariton lasing. These results 
are interpreted as a measurement of the lifetime of the 
Goldstone boson [7J. 

Recent experimental measurements of time-resolved 
PL in similar systems [8] reveal the importance of both 
the dynamics involving a single polarization, and the dy- 
namics involving the two photon polarizations. On the 
other hand, a system with a single photon polarization 
- a single polariton condensate - could be realized by 
means of a strong magnetic field breaking the degener- 
acy between the "right-handed" and "left-handed" con- 
densates. 

In the present paper, I consider a model with a sin- 
gle photon polarization and explore how the way the 



system is pumped and probed influences the decay dy- 
namics of the perturbations. I started from a scheme, 
sketched in Refs. [PJ [TO], in which pumping and photon 
losses in the polariton system are described by means 
of two terms in the master equation for the density ma- 
trix. The linearization of the master equation around 
the stationary solution leads to a system of equations 
with a source term for the small perturbations. The way 
the system is pumped determines the small-oscillation 
modes of these equations, whereas the way the system is 
probed determines which of the eigenmodes are excited. 
In an oversimplified model, with unrealistic parameters, 
I found decay times from 30 to 5000 ps, with no singular 
behaviour at the polariton lasing threshold. The large 
decay times correspond to eigenmodes involving the ex- 
citation of dark polariton states which, under resonant 
pumping, exhibit a bottleneck effect. 

The starting point is a simple expression for the sta- 
tionary spectral function, S(ui), describing the PL emis- 
sion along the cavity symmetry axis, Eq. (19) of Ref. 
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The magnitudes Tu (linewidhts) and luij = (Ej — Ei)/H 
depend only on the many-polariton wavefunctions and 
energies, and the system parameters, P (pumping rate), 
and k = 0.1 ps -1 (photon losses rate). (I\a\J) are the 
matrix elements of the photon annihilation operator be- 
tween the many-polariton wavefunctions |J) and \I). 

In our model, describing a many-exciton quantum dot 
strongly interacting with the lowest photon mode of a mi- 
crocavity |llj . the wavefunctions and energies, and from 
them Tu, lujj, and (I\a\ J), are obtained by numerically 
diagonalizing the electron-hole-photon Hamiltonian. On 
the other hand, the stationary solutions, pf ^ should be 
obtained from the master equation for the occupations 
(coherences are three orders of magnitude smaller |10j 
and will be neglected): 
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N P oi{J) is the polariton number (number of electron-hole 
pairs plus number of photons, which is a conserved quan- 
tity) of state | J), and Pji is the pumping rate from state 
J to state /. Notice that Eqs. Q depends only on k, 
P, and the matrix elements (I\a\J). These equations are 
linearly dependent, which expresses the conservation of 
probability: 



E^ = 1 - 
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The stationary solutions, p\ , are obtained by making 
the l.h.s. of Eqs. ^ equal to zero, and complementing 
this homogeneous linear system with the normalization 
condition, Eq. 

Under a small perturbation of the pumping rate, SP(t) , 
there is a variation of the density matrix, 8p(t), and a 
variation of the spectral function: 
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FIG. 1: Uniform pumping, (a) Position of the PL line as a 
function of the pumping rate, (b) Photon second-order co- 
herence function. 



The response, Spj, to the pulsed perturbation SP(t) 
satisfies the linear system obtained by varying Eqs. 
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which should be complemented with: 
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The structure of Eqs. ([5| is the following: d<5p/dt = 
ASp + 5P(t)Bp(°°> . The eigenvalues of matrix A are the 
small oscillation frequencies of the system, whereas the 
probe perturbation, matrix B, and the stationary den- 
sity matrix conform the source term determining which 
oscillation modes are excited by the probe pulse. 

Let us recall the diagonal terms in Eqs. (|5| in order to 
get a qualitative understanding of the decay times: 



^ = -n^\(J\a\I)\ 2 S Pl 
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This equation shows that the excess occupation of state /, 
dpi, can decay in two ways. The first term corresponds to 
photon emission, and the second to pumping the excess 
occupation to a higher polariton state, J, which may 
further decay through photon emission. When the state 
/ is dark, J^j l(^l a l-0| 2 ~ 0, only the second mechanism 
acts. If, in addition, 8pi is not pumped to higher states 
because of selective pumping (Pjj w 0) , then the decay 
of 5pi may take very long times. Below, we consider 
different regimes of pumping and probing the polariton 
system. 

(a) Uniform pumping and uniform perturbation 

In this case, P u = P, and SP u (t) = 6P(t). This 
situation seems to correspond to laser excitation energies 
well above the lower polariton branch, and perturbations 
at these higher energies. 

I show in Fig. [T] the position of the main PL line as 
a function of the pumping rate, and the corresponding 
photon second-order coherence function, g( 2 >(0). The 
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jump in the position of the line, and the values near one 
of <p 2 '(Q) identify the threshold for polariton lasing at 
P thr w 0.00014 ps -1 in the model, where I use the fol- 
lowing states in order to solve Eqs. ^ for the stationary 
density matrix: the vacuum (/ = 1), the 17 existing one- 
polariton states in the model (1 = 2 — 18), the 256 exist- 
ing two-polariton states (/ = 19 — 274), and 256 states 
in each sector with 2 < N po i < 10. 

I will study the decay dynamics of the probe for P 
values in the vicinity of Pthr- The probe pulse is taken 
in the following way: 



6P(t) = 10 _5 exp-(f-l) 2 ps"\ 
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where t is given in ps. We find the Spi from Eqs. (|5| and 
compute, as in the experiment [5], the energy-integrated 
differential PL response: 



(9) 



The sum over / is restricted to states such that \Ej — 
Ej — E re f\ < SE, where 5E = 2 meV, and the reference 
energy in the present case is E re t = 1489.2 meV. 

I draw in Fig. [2] (a) the computed SS(t) for three val- 
ues of the pumping rate. They show characteristic decay 
times of around 30 ps. These values can be understood 
from the eigenvalues of the linearized decay modes. In- 
deed, let us write the smallest (in absolute value) eigen- 
values in the present case: -0.0358, -0.0287, -0.0117, 
-0.0067, -0.0023 ps" 1 . Notice that they are real, i.e. 
purely decaying (non propagating) modes. The first of 
them (-0.0358 ps -1 , decay time ~ 28 ps, there are many 
such eigenvalues) corresponds basically to an eigenvec- 
tor in which a single dark one-polariton state is excited. 
The decay is due to pumping to two-polariton states with 
further emission of photons. Notice that Pthr = 0.00014 
ps -1 times the number of available two-polariton states 
(256) is equal roughly to the eigenvalue, -0.0358 ps -1 . 
These are the modes dominating the observed behaviour 
of SS(t) in the present case. Let us stress that the eigen- 
values practically do not change when P is varied around 
Pthr- Thus, we can not relate the slowest mode (or any 
other) to the lasing transition. In the scheme of pertur- 
bation I am using, modes with larger decay times are 
not excited. In particular, the last one (-0.0023 ps -1 , de- 
cay time ~ 434 ps) corresponds to an eigenvector involv- 
ing the simultaneous excitation of one- and two-polariton 
dark states combined with higher polariton states. We 
shall see that slower decaying modes can be observed by 
means of a selective excitation. 

(b) Uniform pumping, selective perturbation 
At this point I consider a situation in which the sys- 
tem is pumped at high energies, as above, but resonantly 
perturbed. This means that SPu(t) is given by Eq. dsj , 



o.ooi 



a 0.0008 
xi 

J-H 

3 0.0006 

-a 

73 0.0004 
B 

bO 0.0002 



1 1 ' 1 

A (a) Uniform Pumping, Uniform perturbation 

A 


•.\\ 

" "A 


P=0.00012ps"' 




A 
w\ 


- - P=0.00014ps"' 




"?\ 


- ■ ■• P=0.00015 ps"' 




1,1, 







20 



3e-05 



2e-05 



40 60 

t(ps) 



80 



100 



a 
a 

i 



« -le-05 




If 

2e-05 i 



50 



- - P=0.00014ps" 

- ■ ■• P=0.00015 ps"' 



L 

100 

Hps) 



i 

150 



200 



FIG. 2: (Color online) Time evolution of the energy- 
integrated differential PL response, Eq. Q, under uniform 
pumping, (a) Uniform perturbation, (b) Selective perturba- 
tion. 



only when the energy difference satisfies the inequality 
\Ej — Ej — E pertU rb\ < SE. Otherwise it is zero. We 
show results for E pertur b — 1499.25 meV, and SE = 2 
meV. The chosen E pertur b corresponds to the resonant 
excitation of a dark one-polariton state, labeled by the 
number 15 (E pertU rb = Eis — E\). Of course, other tran- 
sitions may have the same, or close, excitation energy, 
and could be excited if they satisfy the above inequality. 
Below, we shall discuss in more details, how a dark state 
could be perturbed. 

The eigenvalues have not changed because I have not 
modified the pumping scheme. But the energy-integrated 
SS(t) shows decay times of around 60 ps in the present 
case, indicating the excitation of slower-decaying modes, 
as compared to uniform perturbation. Results are drawn 
in Fig. [2] (b), where an oscillation at early times is also 
observed. It is possible that with a different perturbation 
energy, E pertur b, or a different perturbation strategy the 
slowest mode could be reached as well. 

(c) Selective pumping, selective perturbation 
Finally, we consider resonant pumping. Selective 
pumping can increase dramatically the decay times be- 
cause there could be almost dark states, not connected to 
higher states through pumping. The decay of an excess 
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FIG. 3: Selective pumping, (a) Position of the PL line as 
a function of the pumping rate, (b) Photon second-order co- 
herence function. 




1500 
1495 
<U 1490 



I 



£1485 
1480 



1475, 
1470 



(b) 







40 60 

t(ps) 



o-oL = 
mL= 1 



_l I L_ 

10 

State Label 



15 




20 



FIG. 4: (Color online) (a) Time evolution of the energy- 
integrated differential PL response under selective (resonant) 
pumping and selective perturbation. (b) The lowest 2- 
polariton states with L = and L = 1. 



occupation in these states could take very long times. 

The stationary pumping rate is chosen in the form: 
Pu = P only when \Ej - Ej - E pump \ < 8E, where 
Epump — 1488 meV, and 6E — 2 meV. The position of 
the main PL line and the photon second-order coherence 
function are shown in Fig. [3] The jump in the position 
of the PL line and the abrupt variation of <7^(0) allow 
us to identify the threshold rate: P t hr ~ 0.04 ps -1 . No- 
tice that this value is much higher than the threshold 
under uniform pumping, something reasonable. Notice 
also that even below threshold the coherence function 
take values very close to one. This initial coherence is, 
in some sense, inherited from the pumping. 

The energy-integrated PL response is drawn in Fig. [4] 
(a) for P = 0.04 and P = 0.06 ps" 1 . After an initial 
transient period, the curves become almost flat, suggest- 
ing the excitation of very slow decay modes. There are 
many eigenvalues of the linearized problem with absolute 
value lower than 0.01 ps -1 (decay times larger than 100 
ps). The smallest of them is -0.0002 ps -1 , that is a decay 
time of 5000 ps. Let us stress that no singular behaviour 
of the eigenvalues across the threshold is found, which 
means that we can not relate any of the eigenmodes to 
the lasing threshold. 

Let us consider the question about the selective (reso- 



nant) pumping or perturbation of a dark state. Electron- 
hole pairs with a definite energy could be injected to the 
system, but this seems to be difficult to control. On the 
other hand, the direct optical transition is prohibited be- 
cause, by definition, the state is dark. However, notice 
that we are considering transitions between L = states 
[11 , which are responsible for the PL emission along the 
cavity symmetry axis. States with L = 1 or higher, very 
close in energy to L = states, are very common, as 
can be seen, for example, in Fig. Qb). These L = 1 
states could be optically excited (with a non-zero trans- 
ferred linear momentum, k), and then may decay towards 
L = states through emission of very low-energy acous- 
tical phonons. In the reported experiment [S], the sys- 
tem is both pumped and perturbed by using k ^ laser 
beans. 

In conclusion, we studied a model polariton system 
with a single photon polarization and computed decay 
times of probe pulses against a stationary polariton dis- 
tribution. Under non-resonant pumping conditions, the 
computed decay times are 30-60 ps, whereas under res- 
onant pumping very large decay times, of the order of 
thousands of picoseconds, are obtained. The dark polari- 
ton states play a fundamental role in the decay dynam- 



ics, specially under resonant pumping, where the excess 
(probe) occupations of particular dark states may take 
very long times to decay. No singular behaviour of any 
decay mode at the threshold for lasing is observed. 
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